Co-ordination between Rashba spin-orbital interaction and space charge effect and 

enhanced spin injection into semiconductors 
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We consider the effect of the Rashba spin-orbital interaction and space charge in a ferromagnet- 
insulator/semiconductor/insulator-ferromagnet junction where the spin current is severely affected 
by the doping, band structure and charge screening in the semiconductor. In diffusion region, if the 
the resistance of the tunneling barriers is comparable to the semiconductor resistance, the magne- 
toresistance of this junction can be greatly enhanced under appropriate doping by the co-ordination 
between the Rashba effect and screened Coulomb interaction in the nonequilibrium transport pro- 
cesses within Hartree approximation. 
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The spintronics is interesting since this involves explo- 
ration of the extra degrees of freedom provided by elec- 
tron spin, in addition to those due to electron charge, 
which is believed to be very useful in manipulating fu- 
ture electronic devices [1]. To realize such a spin device, 
the Rashba spin-orbit interaction is often considered [2] . 
Since this is caused by structural inversion asymmetry in 
quantum wells, it can be artificially controlled by adjust- 
ing the applied gate voltages and specifically designing 
the hetcrostructure [3]. 

On the other hand, one of the current focus in spin- 
tronics lies in injecting spin polarized electrons into non- 
magnetic semiconductors [4-9]. This is partially moti- 
vated by the high magnetoresistance observed in ferro- 
magnet tunnel junctions [10]. However, in the diffusion 
region and the room temperature, experiments have so 
far observed a small magnetoresistance ratio of 1% [11,12] 
in ferromagnet/semiconductor/ferromagnet structures. 
Rashba proposed a ferromagnetic metal/tunncling- 
insulator/semiconductor (FIS) junction to improve the 
spin injection rate [15]. Many different geometries of the 
tunnehng junction were discussed in Ref. [16]. It is pos- 
sible that to lift the magnetoresistance to near 10% with 
a very high electric field [17]. 

In a recent work [18], we investigated the space charge 
effect in the nonequilibrium transport process within the 
Hatree approximation and found that the magnetore- 
sistance of a double FIS junction could be greatly in- 
creased if one carefully adjusts the parameters of the 
junction, such as the charge screening length and the 
size of the semiconductor. The space charge effect in the 
nonequilibrium transport process first was found play- 
ing an important role in the characteristics of the de- 
vice of ferromagnetic/non-magnetic/ferromagnetic metal 
junction [19]. Under steady state nonequilibrium condi- 
tions, a magnetization dipole layer much larger than the 
charge dipole layer is induced at the interface while the 
magnetization dipole layer is zero under equilibrium con- 
ditions. We have applied this idea to the double FIS 
junction [18]. 

We now ask the question: How do co-ordinations be- 



tween the Rashba spin-orbital interaction and the space 
charge affect the spin transport in the double FIS junc- 
tion? It was known that the Coulomb interaction may 
enhance the Rashba effect. In a one-dimensional Lut- 
tinger liquid formalism, Hausler has reported an enhance- 
ment of the Rashba effect due to the spin charge sepa- 
ration [20]. Such an enhancement was also found in a 
two-dimensional electron gas system [21]. However, the 
effect of the Rashba spin-orbital interaction in the dif- 
fusion region was thought to be not important. In this 
work, we would like to investigate the co-ordination be- 
tween the Rashba spin-orbital interaction and the space 
charge effect in the double FIS junction in the nonequilib- 
rium transport process. To completely solve the nonequi- 
librium problem with interactions is highly non-trivial. 
It requires solving self consistently Boltzmann type spin 
transport equation with the Poisson equation. What we 
want to touch in this work is using a simple Hartree ap- 
proximation to check if the effect comes from these in- 
teractions is significant. Consequently, it is found that 
if (i) the resistance of the tunneling barrier is compa- 
rable to the semiconductor resistance and (ii) the n-type 
semiconductor has an appropriate doping, then while the 
magnetoresistance of the junction is greatly enhanced as 
the charge screening length becomes shorter as we al- 
ready showed in [18], it is grown up as the Rashba spin- 
orbital coupling. The shorter charge screening length 
is , the more gain in the magnetoresistance comes from 
the Rashba term. Comparing with the non-interacting 
model where the Rashba effect may not be observed in 
the diffusion region, we see a great increment of the mag- 
netoresistance comes from the co-ordination between the 
Rashba spin-orbital interaction and space charge effect. 

Rashba Spin- Orbital Coupling: We consider a two- 
dimensional electron gas in a narrow gap quantum well, 
such as those based on InAs. The spin-orbital coupling 
in this kind of systems is dominated by Rashba term [22]: 

Hso = ai<JxPy ~ (^yPx), (1) 

where a^^y are the Pauli matrices and the Rashba param- 
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ctcr a is determined by the asymmetry of the potential 
confining the electrons in the two-dimensional x-y plan 
and can be controlled by a gate voltage [3] . For the sys- 
tem we are considering, its value is about 10~^^ — 10~^^ 
eVm [23]. For a double FIS junction where the current 
flows in the x-direction, we can estimate the spin depen- 
dence of the density of state and the resistance of the 
electron gas. According to the Rashba term (1), the sin- 
gle particle dispersion reads 



2m* 



aakcr , 



(2) 



where m* w 0.04TOe for the bare electron mass rrie and 
ka = k+a^^^; a is called spin-orbital coupling label [24]. 
The Rashba spin-ortial interaction contributes a current 
jso to the system. However, in the diffusion region, such 
an effect may be small [24]. Here, we would like to see 
another effect arising from the Rashba term. When the 
Rashba term is non-zero, it is easy to see the density of 
state NIf is given by [24] 



7V- = 7V^(l + ajv), 



(3) 



where Np is the density of state at the Fermi surface 
of the semiconductor for a = and aj^ ~ for the 
Fermi velocity vp- It is seen that the density of state 
is not cr-dependent [24]. For a typical two-dimensional 
electron gas, the electron density is of order pe = 10^^ 
and 10^^ cm~2. Thus, using vp = ^ = 
is in order of 10~^ and 10~^. In the following, we will see 
that the magnetoresisatnce of the double FIS junction 
is very sensitive to the change of the charge screening 
length which is determined by the density of state. 

Description of the Junction: In the diffusion region and 
the room temperature, we consider the junction Fl-Il-S- 
I2-F2 along the .x-direction. The thickness of the metal 
(Fl, F2), the insulating barriers (II, 12) and the semicon- 
ductor (S) are denoted by L^'^, di^2 and a;o, respectively 
(See, Figure 1). For the practical case, xq is less than the 
spin diffusion length I n of the two-dimensional electron 
gas in the semiconductor. The charge screening lengths 
in the metal and the semiconductor are denoted by Al^^ 
and Xn, respectively. In our model, we assume Aat <C /jv 
and Xr^l <C Ir^l, the spin diffusion lengths in the metal. 
Typically, A ~ 10~^nm, / ~ lO^nm, and In > 1/xm. 
xq ~ lOOnm to l/xm, depending on the structure of junc- 
tions. The screening length in the semiconductor, A^v, is 
dependent on the doping of the semiconductor and can 
vary in a wide range, say lOnm for the heavy doped semi- 
conductor and lOOnm-l/im for the lightly doped or un- 
doped semiconductor. To avoid discussing the spin-orbit 
splitting of the heavy and the light-holed bands near the 
zone center, we consider the n-type semiconductor only. 

The problem we would like to solve has been defined 
in [19,18]. Summarily, it is four sets of equations 

(1) The total charge-current conservation is described by 



(4) 



where p is the charge density. 

(2) The spin s-dependent current is determined by the 

diffusion equations. In order to allow an analytic analy- 
sis, we take a simple Hartree approximation into account 
to see the space charge effect. Under such an approxima- 
tion, the diffusion equations reads [19] 



(5) 



where the magnitude of the electric charge has been 
set as one: E is the external electric field; the chemi- 
cal potential /i., is related to the charge density ps by 
V/is = where Ns is the spin-dependent density of 
state. Wo = J dfUintir — f')p{f') is the potential caused 
by the screened Coulomb interaction Uint{r) and VVFq is 
the so-called screening field induced by Uint{r)- 
(3) The third is the magnetization relaxation equation 
where the magnetization density M = p^ — pi relaxes 
with a renormalized spin diffusion length /. In the relax- 
ation time approximation, one has 



V^M - M/l^ = 0. 



(6) 



(4) The boundary conditions at the interfaces are given 

by 



Afis - AW = r(l - sj)js 



(7) 



where ps = + Ex; Ex is the voltage drop on the left 
side of the barrier and AW is the electric potential drop 
across the barrier, which is assumed much smaller than 
A/is; = ^(1 — S7) is the barrier resistance. We assume 
that there is no spin relaxation in the insulator, the spin- 
dependent currents are continuous across the junctions, 
j^(-di/2) = if (di/2) and j^{xo + (1^/2) = j^{xo - 
d2/2). 

In addition, we have the neutrality condition for the 

total charges {Ql.n.r, e.g., Qm = l^l^-i^^"^ P'^^) accu- 
mulated at the interfaces. By Gauss' law, for the point 
d\l2 < X < xo + di/2, i.e., inside the semiconductor, the 
potential Wq is determined by 



f 

Jdi 



/2 



pdx, 



(8) 



whose constant part of the right hand side gives the con- 
straint on the charge while the x-dependent part gives the 
function form of the potential Wq. Another constraint is 
the neutrality of the system: 



Ql + Qn + Qr = 0. 



(9) 



With those sets of equations (Eqs.(4)-(7)) and the two 
constraint on the charges (Eqs.(8) and (9)), the problem 
can be solved. The formal solutions of the problem are 
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''L 



(10) 



with similar solutions for the right hand side. In the 
semiconductor, if Ajv <C a;o, 

p(2)(^) = ^p(2)e(--o+c^./2)/A„ 

In 



X 



i^/,(^)p(^-^o+'^2/2)/;iv 



+ 72 P20 ^ 
'at 



(11) 



can be obtained similarly. All of coefficients in 
(10) and (11) can be determined by using Eqs.(4)-(7) 
and the constraint (8) and (9). The screening potential 
Wo is determined by Gauss' law. The total current is 
j = J2sjs- Although js are not a constant, the total 
current j is still a constant. 

The spin- dependent currents: It is necessary to simplify 
the problem to demonstrate the essential physics. One 
sets the parameters of the metals and barrier widths on 
the left and right sides to be the same: Xj^ = Xl = X, 
II = Ir = I, d\ = d2 = d and so on. The resistances of 
the barrier layers are taken as r^^^ = r^^^ = t; 71 = 72 = 
7 for the parallel configuration and 71 = —72 = 7 for 
the anti-parallel configuration. To illustrate, we focus on 
the calculation in the left barrier located at a; = 0. The 
tunneling resistance is given by 

r^P = r(l - 7s) = ro,s exp[d{Ks{p) - Ks{0))], (12) 

where Ks{p) oc dx[2m{U — Aps{0)x/d)]^/'^ , with U the 
barrier height. The current j^{x) at .t = is dependent 
on the bias voltage and the interaction, which is given by 

j^{0)=AJos, (13) 
where, for the ferromagnet metal on left hand side, 

/ 1 — ap ti 

and jos = (JsE is the current with no interaction and the 
current jso that is contributed to from the Rashba term 
has been neglected; (i {a) measures the spin asymmetry 
of the conductivities Us (densities of states at the Fermi 
surface Ns): = f (1 + /3s) and = ^Nf{1 + sa) 
where Np is related to the screening length A by -j^ = 

27rA^ lZaf3 ■ Noting that both pfg and Mq are propor- 
tional to the external electric field E, Ag is solely deter- 
mined by the material parameters. 



Eq. (13) implies that the spin-dependent current js{0) 
passing the interface differs a factor from the non- 
interacting current jos(O). For the parallel configuration, 
jos(O) is given by [18] 



Jo.(0) 



V 



R'j^ + 2roM0)' 



(15) 



where y,(0) = e''o,d[^^ii-ii-Ap,io)r'^)-i] ^so ^ 

iijv(l — Pn) with /3jv ~ cun- For the anti-parallel config- 
uration, 



V 



(16) 



The magnetoresistance: Since the total current is con- 
stant everywhere, we have, for the parallel configuration 

= jsiO)/V and for the anti-parallel configuration, 

flX^ ~ Ssis(O)/^- From these, we obtain the magne- 
toresistance ratio 



Rap — Rp 



X 



(Rap + Rp) 2 + 



(17) 



1. For non- 



^ — A^ = 1 and we knew that 



where X = 

interacting electrons, A^ 
that AR/R will not be beyond a maximal value (at 
V -^0) about 3.2% for ro+ : ro- : i?jv = 1 : 2 : 1 and de- 
cays as the bias voltage increases [18]. Since jso has been 
neglected, there is no observable Rashba effect without 
interactions. After the interaction is included, the ratio 
(17) grows greatly and the Rashba eff'ect is enhanced as 
Xn becomes smaller. This can be clearly seen in Figs. 2 
and 3. We take xq — 1.25/im, A = O.lnm, I = 20nm, and 
In = 3/xm and set a — P ~ 1/2. The different choice of a 
and /3 will not qualitatively affect the result if they do not 
deviate from 1/2 too much. In Fig. 2, we depict the mag- 
netoresistance versus the bias voltage for Aat lOOnm 
with ajv from to 0.05. It is seen that AR/R raises 
from 16% to 18% when ajv from to 0.05. In Fig. 3, 
for Aat = 50nm, it is shown that the Rashba effect lifts 
much faster. AR/R raises from 25% to 30% for ajv from 
to 0.05. Hence, we see a strong co-ordination between 
the Rashba spin-orbital and the screened Coulomb in- 
teractions to increase the magnetoresistance. While the 
Coulomb interaction largely enhances the magnetoresis- 
tance, the Rashba spin-orbital interaction may enhance 
AR/R. This kind of Rashba effect becomes more signif- 
icant in a shorter charge screening length. 

In conclusions, we have shown the co-ordination be- 
tween the Rashba spin-orbital and screened Coulomb 
interactions on electron injection from ferromagnet to 
semiconductor. The magnetoresistance grows fast as the 
charge screened length in the semiconductor becomes 
shorter. The Rashba term also enhances the magnetore- 
sistance and plays more important role as the Coulomb 
interaction is stronger. In fact, by a chose examination 
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to our solution, this Rashba effect is corresponding to 
a renormalization to the screening length Ajv ^ = 
Ajv/(1 + cxn) in the dominate terms of the solution. Be- 
cause our solution is very sensitive to Xn, it is understood 
that why a small spin-orbital coupling can cause a rela- 
tive large gain in the magnetoresistance. 

This work was partially supported by the NSF of 
China. 
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Fig. 1 The sketch of the double FIS junction. 

Fig. 2 The magnetoresistance versus the bias voltage (in 
unit U) with Aat = lOOnm. The junction parameters arc 
Xq = 1.25/Ltm, A = O.lnm, I = 20nm, and In = 3/xm. 
Three different values of ajv are taken. 

Fig. 3 The magnetoresistance versus the bias voltage 
with \n = 50nm. Other parameters are the same as 
those in Fig. 2. 
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